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I present the factorization(s) of the Wheeler-DeWitt equation for vacuum FRW minisu- 
perspace universes of arbitrary Hartle-Hawking factor ordering, including the so-called 
strictly isospectral supersymmetric method. By the latter means, one can introduce an 
infinite class of singular FRW minisuperspace wavefunctions characterized by a Darboux 
parameter that mathematically speaking is a Riccati integration constant, while physi- 
cally determines the position of these strictly isospectral singularities on the Misner time 
axis. 

In cosmology one can encounter arbitrary parameters that are difficult to de- 
termine by experiment and theory. One famous example is Einstein's cosmological 
constant, which does set the time scale in the vacuum dominated universe. 1 An- 
other, very recent example is Immirzi's parameter in quantum general relativity. 2 
In this work, using the factorization(s) of the Whcclcr-DcWitt (WDW) minisuper- 
space equation for the vacuum Friedmann-Robertson- Walker (FRW) universes, an 
infinite class of singular FRW wavefunctions are introduced containing an arbitrary 
parameter of Darboux type 3 determining the position of the singularities during 
their Misner time evolution. 

Recalling that factorizations of second order linear one-dimensional differential 
operators are common tools in Witten's supersymmetric quantum mechanics 4 and 
imply particular solutions of Riccati equations known as superpotentials, I explore 
here the physical consequences of the most general factoring of the FRW WDW 
equation for arbitrary Hartle-Hawking 5 factor ordering Q by means of the general 
solution of Riccati equation, a procedure which has been first used in physics by 
Mielnik for the quantum harmonic oscillator, 6 and which is known as the strictly 
isospectral supersymmetric method and/or the double Darboux method. 7 Stated 
differently, I shall exploit the non-uniqueness of the factorization of second-order 
linear differential operators on the simple example of the FRW WDW differential 
equation with Hartle-Hawking factor ordering, that I write down in the form 
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The independent variable fl is Misner's time parameter related to the volume of 
space V at a given cosmological epoch through Q = — lnfV 1 / 3 ). 8 Using the change 



of function ^ = <&exp(— ^Q), one gets the following Schrodinger equation of zero 
eigenvalue on the Misner axis 

$" - (91 + e -4nA $ = o (2) 



The general solution is the modified Bessel function Zq/ 4 (z ie~ 2S1 ) = CilQ/4(^e~ 2n )^ 
C2KQ/ i (^e~ 2Q ) (the Cs are superposition constants). Usually, the strictly isospec- 
tral supersymmetric technique requires nodeless, normalizable Schrodinger solu- 
tions. Although the modified Bessel functions Zq/4 are not normalizable, one can 
still use the strictly isospectral technique as first shown by Pappademos et al for a 
few simple scattering examples in ordinary quantum mechanics. 9 Using the modi- 
fied Bessel function Zq/^, the strictly isospectral supersymmetric construction says 
that the class of strictly isospectral potentials for the FRW cosmological models are 
given by 

s iso (n-,n D ) = s{n)-2[in(j z {n) + n D )]" 

where S(£l) = ^- +e~ 4 °, fl^ is the Darboux parameter introduced by this method 
and 

Jz(Q) = I Z 2 Q/ Se- 2y mdy . (4) 

J oo 

To get (3) one factorizes the one-dimensional FRW Schrodinger equation with the 
operators A = + W(fi) and A* = -J^ + W(f2), where the superpotential 
function is given by W = —Zq^/Zq/^. The Schrodinger potential and Witten's 

superpotential enter an initial 'bosonic' Riccati equation S = W 2 — W . On the other 
hand, one can build a 'fermionic' Riccati equation S + = W 2 + W , corresponding 
to a 'fermionic' Schrodinger equation for which the operators A and A + are applied 
in reversed order. Thus, the 'fermionic' FRW Schrodinger potential is found to be 

S + = ^- +e~ 4 ° — 2( Z Q/ * ) . This potential does not have Zq/^ as an eigenfunction. 
However, it is possible to reintroduce the -Zq/4 solution into the spectrum, by using 
the general superpotential solution of the fermionic Riccati equation. The general 
Riccati solution reads 



w gen = w(o) + ^ in[j z (fi) + n D ], i .-, ) 



where Hp € (— oo, +oo) occurs as an integration constant. 10 The way to obtain (5) 
is well known since the work of Miclnik and will not be repeated here. From it one 
can easily get (3). The Darboux parameter family of FRW wavefunctions differs 
from the initial one, being the ^-dependent quotient 11 

Jz + ^d 



z gen {n- n D ) = T\ 2 . (6) 
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The last equation has interesting implications for the celebrated issue of cosmolog- 
ical singularities. Indeed, the strictly isospectral supersymmetric method may be 
considered as a way of generating parametric families of new cosmological wavefunc- 
tions with the important feature that there is always a point Q s on the Misner axis 
defined by flo = —Jzi^s) at which they become singular (as well as the strictly 
isospectral potentials). In other words, an infinite set of such strictly isospectral su- 
persymmetric singularities occur, whereas in the classical case one is usually bound 
only to the idea of the original singularity, which is by definition at Q — > oo, "an 
infinite fi time in the past" as Misner put it. 8 Therefore, we agree with Wheeler's 
dictum 12 "No, quantum mechanics does not provide an escape from gravitational 
collapse" . Even more, the solutions with singularities are never a set of measure 
zero for the quantum FRW universes. There exist well-defined mathematical pro- 
cedures to produce an infinite amount of solutions with singularities, at least for 
the FRW case. Finally, for some plots displaying the strictly isospectral supersym- 
metric FRW singularities in the particular case Q = the reader is directed to one 
of my works with Socorro. 13 
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